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LIST OF CURRENT CLAIMS 

1. (currently amended) A statistical facial feature extraction method, comprising: 
a first procedure for creating a statistical face shape model based on a plurality of 
training face images, including: 

an image selecting step, to select N training face images; 

a feature labeling step, to respectively label feature points located in n different 
blocks of the training face images to define corresponding shape vectors of the training 
face images; 

an aligning step, to align each shape vector with a reference shape vector to thus 
obtain aligned shape vectors; and 

a statistical face shape model computing step, to use a principal component 
analysis (PCA) process to compute a plurality of principal components based on the 
aligned shape vectors to create a statistical face shape model, wherein the statistical face 
shape model represents the shape vectors by combining a plurality of projection 
coefficients , and the statistical face shape model computing step includes: computing a 
mean value of the feature points of the aligned shape vectors to define a mean shape 
vector x , subtracting each aligned shape vector x n by the mean shape vector x _ to form a 
matrix A, computing a covariance matrix C of the matrix A, and computing a plurality of 
principal components according to eigenvectors to form the statistical face shape model ; 
and 

a second procedure for extracting a plurality of facial features from a test face 
image, including: 

a test face image selecting step, to select a test face image; 

an initial guessing step, to guess initial positions of n test feature points located in 
the test face image, wherein the initial position of each test feature point is a mean value of 
the feature points of the aligned shape vectors; 
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a search range defining step, to define n search ranges in the test face image, based 
on the initial position of each test feature point, wherein each search range corresponds to 
a different block; 

a candidate feature point labeling step, to label a plurality of candidate feature 
points for each search range; 

a test shape vector forming step, to do combination of the candidate feature points 
in different search ranges in order to form a plurality of test shape vectors; and 

a determining step, to match the test shape vectors respectively to both the mean 
value and the principal components for computing a similarity, and to accordingly assign 
one feature point corresponding to one, having the best similarity, of the test shape vectors 
as facial features of the test face image. 

2. (original) The method as claimed in claim 1, wherein in the feature labeling step 
of the first procedure, the feature points are coordinates for corners of eyes and mouth on 
each training face image. 

3. (original) The method as claimed in claim 1, wherein the feature labeling step of 
the first procedure manually labels the feature points of each training face image. 

4. (original) The method as claimed in claim 1 , wherein the reference shape vector 
is one of the shape vectors. 

5. (original) The method as claimed in claim 1, wherein the aligning step of the 
first procedure uses a 2D scaled rigid transform algorithm to align each shape vector with 
the reference shape vector. 



6. (original) The method as claimed in claim 5, wherein the aligning step of the 
first procedure further comprises the steps of: 
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selecting the reference shape vector as x. = (x iI ,y iI ,...,x ill ,y ill ) and one of the 
shape vectors asx^ (x , y., x jn , y jn ) ; 

computing a squared Euclidean distance E between the vectors Xj and Xj based on 
the following equation E = (x ; — M (N) (a,0)[x j ] — t) T (X; — M (N) (a,0)[x j ] — t), 

where M (N) (a,0)[x j ] - t is a geometric transform function defining with a plurality of 
transfer parameters to align the vector xj, the transfer parameters include a rotating angle 
0 , a scaling factor a , and a shifting vector represented by t = (t x ,t y ) , and M (N) (a,0) 
is a 2nx2n diagonal blocked matrix as well as 
( acos0x jk -asin0y jk ^ 



M(a,0) 



|_y jk J ^asin0x jk +acos0y jk 
fa cos 0 -a sin Q\ 

M < a > 6 H • A A F 



for 1= k= n, as 



finding the smallest squared Euclidean distance and corresponding rotating angle 
0 j5 scaling factor a } and shifting vector represented by tj = (t xj , t yj ) to align the shape 

vector Xj so similar as the reference shape vector x<; 

computing a sum of smallest squared Euclidean distances after the N shape vectors 
are all aligned so similar as the reference shape vector, ending the aligning step when the 
sum is smaller than a predetermined threshold; 

computing a mean value of the feature points in each block for the aligned shape 

1 n 

vectors to define a mean shape vector for each aligned shape vector as x = — £ x a , 

wherein x a is the aligned shape vector; and 

assigning the mean shape vector as the reference shape vector and the aligned 
shape vectors as the shape vectors and then repeating the aligning step until all shape 
vectors are aligned. 
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7. (original) The method as claimed in claim 6, wherein the transfer parameters is 
obtained by a least square algorithm. 

8. (original) The method as claimed in claim 1, wherein the statistical face shape 
model is a point distribution model (PDM). 

9. (currently amended) The method as claimed in claim [[8]]I, wherein in the 
statistical face shape model computing step, the mean shape vector x _ is defined as 
- 1 * 

x = — IX 

j_j a ~' , where xa is the aligned shape vector; the matrix A is formed as 
A = |d xi ,d X2 ,...,d XN J J _where_d x = x a — x ; the covariance matrix C is formed as 

C = AA T ; the eigenvectors are derived from Cv° = A? k V k with eigenvalues 

corresponding to the covariance matrix C, where Af,, represents eigenvalues of the 

covariance matrix C, V* represents eigenvectors of the covariance matrix C, and 1= k= 

m, m being the dimension of the covariance matrix C for A, 5 , >X 2 > ... > X m th e statistical 

fac e shap e mod e l computing step of th e first procedure furth e r compris e s th e st e ps of: 

computing a m e an valu e of th e f e atur e points of th e align e d shap e v e ctors to d e fine 
- 1 N 

a m e an shap e v e ctor as ^ a_1 , wh e r e in xa is th e align e d shap e v e ctor; 

subtracting e ach align e d shap e v e ctor by th e m e an shap e v e ctor to form a matrix as 

A ~ K A 2 ? -» d x N J» wherein d Xa = x, — *i 

computing a covarianc e matrix of th e matrix A as C = AA 1 ; and 
computing a plurality of principal compon e nts according to e ig e nv e ctors deriv e d 
from Cv* — A, k v k with e ig e nvalu e s corr e sponding to th e covariance matrix C, to form 
th e statistical fac e shap e mod e l, wh e r e in - ?c^- r e pr e s e nts e ig e nvalu e s of th e covarianc e 
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lv e otors of the 



matrix C, and 1 ^ k 4 m, wher e m 



th e dim e nsion of th e 



matrix C for X s , > A, s 2 > ... > A, s m . 



10. (original) The method as claimed in claim 1, wherein each shape vector xj 
consists of n feature vectors Sj located in different blocks, so an average value as 

1 n 

t = — X s of the feature vectors Sj corresponding to special blocks of all shape vectors xj 
Nh j 

is defined as a feature template. 

11. (original) The method as claimed in claim 1, wherein in the initial guessing 
step of the second procedure, scaling of initial guess shapes formed by the test feature 
points is aligned so similar as the test face image. 

12. (original) The method as claimed in claim 10, wherein the candidate feature 
point labeling step of the second procedure further comprises the steps of: 

k 

labeling a plurality of reference points derived from I. = t + Y d b i p j respectively 

j=i 

in each search range, where t is the feature template of block corresponding to a search 
range, pj is j-th principal component of the statistical face shape model computed from the 
training feature vectors, and bj is an associated projection coefficient; 

k 

using £ =|| I ; — t — Xb p • || 2 to compute an error value between one of the 

j=i J 

reference points and the corresponding principal component pj and projection coefficient 
bj; and 

selecting preceding k smallest error values and defining the k smallest error values 
as feature points of the search range. 
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13. (original) The method as claimed in claim 12, wherein the test shape vector 
forming step of the second procedure does combination of the candidate feature points in 
different search ranges to thus form k" test shape vectors. 

14. (original) The method as claimed in claim 10, wherein the determining step of 
the second procedure further comprises the steps of: 

using the average value of the aligned shape vectors and the principal components 
of the statistical face shape model to represent an approximate value of the test shape 

vector as x = x + £bj X p* , where X is a mean shape vector defined according to the 

j=i 3 

mean value of the feature points of the aligned shape vectors, p* is j-th principal 
component of the statistical face shape model, and bj is a corresponding projection 
coefficient; 

using a 2D scaled rigid transform algorithm to align the test shape vector 

represented by X = M(a,0)^x + Zb* p*J + t , where 0 , a and t are a rotating angle, a 

scaling factor and a shifting vector respectively; 

computing a normalized distance of the aligned test shape vectors by 



assigning one candidate feature point corresponding to one, having the smallest 
normalized distance, of the aligned test shape vectors as facial feature of the test face 
image. 

15. (original) The method as claimed in claim 10, wherein the determining step of 
the second procedure further comprises the steps of: 

using the average value of the aligned shape vectors and the principal components 
of the statistical face shape model to represent an approximate value of the test shape 
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vector as x = X + Zbjp- , where X is a mean shape vector defined according to the 

mean value of the feature points of the aligned shape vectors, p* is j-th principal 

component of the statistical face shape model, and bj is a corresponding projection 
coefficient; 

using a 2D scaled rigid transform algorithm to align the test shape vector 

represented by x = M(a,0)j^x + Zbj P*J + t , where 0 , a and t are a rotating angle, a 

scaling factor and a shifting vector respectively; 

computing an error value between the test shape vector and the average value of 

n k 

the aligned test shape vectors by e(x) = w, £|| I. (x) - tj - ZbjP- 1| 2 +w 2 d(x) , where 

i=i j=i 

n k . . 

XII I (x) — tj — SbjP || 2 i sa similarity of the test shape vector to corresponding aligned 

i=i j=i 1 

shape vector x a , and d(x) is the normalized distance of the aligned test shape vectors; and 

assigning one candidate feature point corresponding to one, having the smallest 
error value, of the test shape vectors as facial feature of the test face image. 



16. (original) The method as claimed in claim 15, wherein the error value is 

f- fWV 

computed by an equation e(x) = W |]£ — + W d(x) . 



